
function z=MOP2(R)

D=5; A=D/R(2); C= D/R(3); B= sqrt(A^2+C^2+D^2-2*R(1)*A*C);
% r=0;
%     z1=abs(((1/x(3))-(x(1)/x(2))+(x(3)/x(2)))/sqrt(1+(1/x(2)^2)+(1/x(3)^2)-(2*x(1)/(x(2)*x(3)))))-sin(r);
%     z2=abs(((1/x(3))-(x(1)/x(2))-(x(3)/x(2)))/sqrt(1+(1/x(2)^2)+(1/x(3)^2)-(2*x(1)/(x(2)*x(3)))))-sin(r);
z1=abs(pi/2 - acos((C.^2 + B.^2 - A.^2 - D.^2 + 2*A*D)./(2*B*C)));
z2=abs(acos((C.^2 + B.^2 - A.^2 - D.^2 - 2*A*D)./(2*B*C))- pi/2);
z=[z1;z2];
    
end
